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FOREWORD 

This  report  describes  a  method  for  determining  the  fillet  radius 
likely  to  prevent  separation  of  the  corner  flow  at  the  intersection  of 
an  airfoil  and  a  hub  or  shroud  where  no  relative  motion  is  present. 

The  work  was  performed  in  the  Turbine  Engine  Division  of  the 
Aero  Propulsion  Laboratory,  Air  Force  Wright  Aeronautical  Laboratories, 
Wright-Patterson  Air  Force  Base,  Ohio.  It  was  accomplished  under  Project 
2307,  Task  SI,  Work  Unit  27.  The  effort  was  conducted  by  Dr.  Lucien  L. 
Debruge  (AFAPL/POTX)  during  the  period  April  1980  to  February  1981. 
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SECTION  I 
INTRODUCTION 

Based  on  an  analysis  presented  in  Reference  I,  this  report  describes 
a  procedure  to  calculate  the  parameters  which  determine  the  inception  of 
separation  at  the  intersection  of  a  turbocompressor  stator  blade  row 
airfoil  and  its  adjacent  platform  or  endwall.  These  parameters  are  then 
used  to  determine  the  geometry  of  a  fillet  which  will  prevent  corner  flow 
separation  from  occurring. 

The  report  includes  a  computer  program  which  was  developed  to 
calculate  these  parameters  from  a  knowledge  of  blade  row  inlet  and  exit 
flow  conditions  and  geometry. 

An  example  of  a  typical  calculation  procedure  of  a  corner  fillet 
geometry  is  presented. 
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SECTION  II 
THEORY  AND  METHODS 


1 .  FLOW  MODEL 

Under  consideration  is  the  flow  development  in  a  compressor  stator 
blade  row  in  the  vicintiy  of  the  corners  formed  by  the  blade  suction  and 
pressure  sides  and  endwalls.  The  flow  is  assumed  turbulent  through  the 
blade  row. 


The  inviscid  flow  velocity  in  the  channel  may  be  calculated  either 
from  the  method  developed  in  Reference  2,  assuming  a  circular  arc  blade 
(Figure  1),  or  directly  obtained  from  a  polynomial  interpolation  of  data 
points  (Figure  2). 


Under  the  assumption  of  a  circular  arc  blade,  the  free  stream 
velocity  is  expressed  as 

Ug  =  tj>l  (a)  +  \p2  (a)  n» 

where  Ue  is  the  velocity  along  a  streamline,  n  is  in  the  Y  direction,  and 
Ugx  and  Ugy  are  the  X  and  Y  components  of  Ug  (Figure  3). 

Furthermore, 

\p2  (a)  =  -it  (AUgy/2ST)  sin  (ira/ST), 

where  AUgy  =  Ug^  -  Ug^  and  U^  and  Ug^  are  tangential  components  of 

the  velocity  at  the  entrance  and  exit  of  the  channel,  respectively.  ST 
is  the  length  of  the  flow  path  in  the  blade  row  and  n  =  0  defines  the 
midchannel  streamline  so  that  the  suction  and  pressure  corners  corre¬ 
spond  to  -gs  and  +np,  respectively. 
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With  the  subscript  o  designating  mid-channel  streamline  conditions, 
we  have: 

U  =  constant  =  U  sin  d>, 
exo  eo  Y 

where  <t>  =  arctg  (U  /U  ).  ■  4 

T  3  ex  ey' 

Also,  \ 

Ueyo  -£  <°>  do  +  Ueyol  '  Ueyol  '  (iV2)  0  '  C*S  (to/ST))’ 

Thus , 

v2 

U  =  (U  2  +  U  2  )  '  =  Ip,  (or) 
eo  v  exo  eyo'  ' 


and 


ttAU. 


u  =  u 


eo 


2ST 


sin 


TTO 

ST  n' 


(1) 


The  expression  for  the  velocity  gradient  is,  therefore. 


3U, 

8a 


AU  r2*  sin<»  si"  (f)  jUeyl  + 

/  V  t 


All. 


ey 


ceps 


•(«)'  «(») 


(S)  -]  J 


(2) 


If  interpolation  of  velocity  data  points  is  chosen,  the  expression  for 
the  free  stream  velocity  and  its  first  derivative  in  the  channel  is: 


u  =  Cl  +  C2  X  +  C3  X2  +  C4  X3  +  C5  x" 
ue 


(1 1 ) 


and 


—  =  C2  +  2C3  X  +  3C3  X2  +  4C5  X3. 
dX 


(2') 
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2.  CALCULATION  OF  TWO-DIMENSIONAL  BOUNDARY  LAYER  MOMENTUM  THICKNESS 


The  two-dimensional  flow  boundary  layer  momentum  thickness  used  in 
the  evaluation  of  the  shape  factor  HI 2  and,  therefore,  of  $  and  n>  is 
obtained  from  Truckenbrodt' s  quadrature  (3).  The  expression  for  the 
two-dimensional  boundary  layer  momentum  thickness,  62 >  is: 


<5  2  (a) 


3+2 

TT 

da 


n _ 

n+l 


(3) 


In  this  equation,  n  =  6. 


The  shape  factor,  H12,  is  obtained  from  Garner's  equation  (3): 


5(H12-1.4) 

-e 


62  £Ue 
ue  da 


+  .0135  (H12-1.4) 


(4) 


The  boundary  layer  velocity  profile  power  law  exponent,  n,  is  directly 
obtained  from 


H1 2  = 

n 

With  the  assumption  (1)  that  the  ratio  of  the  three-dimensional  corner 
boundary  layer  extent  a  to  the  two-dimensional  boundary  layer  thickness 
eo  is  approximately  equal  to  2,  B  is  obtained  from 


The  friction  coefficient,  Cf,  is  obtained  from 
Cf  =  0.013  /  (Ue62/v)  ^ 
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3.  CALCULATION  OF  CORNER  FLOW  PARAMETERS 

The  fact  that  secondary  flows  are  always  generated  in  the  endwall 
boundary  layer  precludes  using  the  two-dimensional  boundary  layer  thick¬ 
ness  directly  derivable  from  Equation  3  as  a  basis  for  the  calculation 
of  corner  flow  parameters.  Also,  the  presence  of  secondary  flows  causes 
the  corner  boundary  layer  to  be  markedly  unsymmetrical .  Nonetheless, 
by  assuming,  in  the  vicinity  of  the  corner,  the  boundary  layer  thickness 
on  the  blade  pressure  and  suction  sides  to  be  equal  to  that  resulting 
on  the  endwall  from  the  presence  of  secondary  flows,  the  three-dimen¬ 
sional  corner  boundary  layer  becomes  symmetrical.  The  calculation  of 
the  interference  displacement  thickness  in  the  corner  conducted  according 
to  the  procedure  developed  in  Reference  1,  yields  a  value  of  this  para¬ 
meter  which  may  be  considered  as  an  upper  limit.  According  to  Reference 
1,  this  upper  limit  will  lead  to  a  conservative  value  of  the  fillet 
radius  likely  to  prevent  separation  in  the  corner. 


The  qualitative  nature  of  the  results  sought  confers  a  particular 
interest  on  the  method  used  by  P.  Suter  to  calculate  the  corner  momentum 
thickness  in  Reference  2.  Figure  3  shows  that  the  development  of  the 
corner  momentum  thickness  along  the  blade  pressure  and  suction  sides  is 
virtually  independent  of  the  presence  of  secondary  flows  generated  in 
the  channel.  This  remark  affords  a  significant  simplification  in  the 
calculation  of  the  interference  displacement  thickness  since  the  para¬ 
meters  of  the  turbulent  endwall  boundary  layer  with  secondary  flow  in 
the  vicinity  of  the  corner  may  be  evaluated  from  a  linear,  first  order 
differential  equation, 

a<59„ 

TT  M1  +  Szc  C2  +  M2)  '  Cf>  (6 


instead  of  a  system  of  two  non-linear  partial  differential  equations, 
Reference  2 


Iy  <Ue£  *2)  + 


3Y  [Uex  Uey  ^ 2  +  Dh^1  +  H12  6  ( U  5Uex  +  «  9Uex  \ 

L  J  2  \  ex  3X  ey  9y  / 


3Uex  Cf 


B11  h  Uey  3  y  2  U<2X  Ue 


(7) 


5 


AFWAL-TR-81 -2020 


9_ 

3X 


“exUey  <S2  '  B21*> 

+  Iy 

Uey  <«2  -  Ah  -  fi) 

_  — 

2 

6  H12  (U  +  Uey  -  Bh  UGV  ^ 

2  ex  gx  y  gy  1  1  Gy  g  Y 


?  Ue^ey 


1  + 


m 


(2m+l) (m+l)  <$ 


2j 


(7) 


In  solving  Equation  6,  the  parameters  Ml  and  M2,  as  well  as  the  friction 
coefficient  are  considered  to  vary  along  the  streamline. 

The  three-dimensional  corner  boundary  layer  parameter  62c  obtained 
from  Equation  6  may  be  defined  as  the  average  corner  momentum  thickness 
at  a  particular  location  along  the  corner.  The  definitions  and  hypotheses 
used  in  the  derivation  of  Equation  6  by  P.  Suter  directly  yield  the  result 
that  Xc,  the  coordinate  of  the  intersection  of  the  corner  bisector  with 
the  corner  boundary  layer  surface,  may  be  expressed  as 

Xc  =  S2c/M3. 


The  thickness  and  the  momentum  thickness  of  the  endwall  boundary  layer 
with  secondary  flow,  in  the  corner,  can,  in  turn  (From  Reference  1)  be 
written: 


eo 


and 

(8) 

6  =  eo  (H12  -  1)/(H12  (H12  +  l)). 

2 

According  to  Reference  1,  the  interference  displacement  thickness,  63  ,  is: 


I 


I 

-  -  -(l 
(e-l)(n+ip 


2n2  3 


Xc 


Xc  “  (n+l) [p(2n+l)+l] 


[B(2n+l)+l] 

nS  + 


3-1 


-  Xc 


\ 


d 


_  2eo3 

2n+l  ”  n+l 


% 


(9) 
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4.  DETERMINATION  OF  FILLET  RADIUS  LIKELY  TO  PREVENT  CORNER  FLOW 
SEPARATION 

As  explained  in  Reference  1,  the  maximum  value  of  5  along  the 
suction  or  pressure  corner  is  used  directly  to  determine  the  fillet 
radius  likely  to  prevent  corner  flow  separation.  For  this  purpose, 
Figure  5  reproduced  from  Reference  1  is  used.  The  procedure  may  be 
described: 


(a)  For  the  Reynolds  number  applicable  to  the  channel  flow  under 

consideration,  obtain  directly  from  Figure  4  a  "safe"  value  of  the 

6 

maximum  corner  interference  displacement  thickness.  For  Re  =  10  ,  for 

instance,  <S  ..  could  be  chosen  between  .007  and  .01. 

3  /  L 

(b)  From  the  maximum  value  of  calculated  from  Equation  9  and 
the  "safe"  value  of  <$3  chosen  in  (a),  estimate  the  percent  decrease  in 
6  which  will  prevent  corner  flow  separation. 

(c)  Use  the  desired  percent  decrease  obtained  in  (b)  to  estimate 
from  Figure  5  the  desirable  fillet  radius  to  be  provided  in  the  corner 

from  the  channel  entrance  to  the  location  of  maximum  5  . 

3 


5.  REMARKS  CONCERNING  THE  CALCULATION  OF  H12y  e,  AND  EU 


The  average  value  of  the  Euler  number, 


is  calculated  along  the  mid-channel  streamline,  XMN,  the  upper  limit  of 
the  integral,  is  the  value  of  X  at  which  the  minimum  value  of  the  mid¬ 
channel  free  stream  velocity  occurs. 
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On  an  airfoil  or  in  a  constant  pressure  gradient  diffuser,  entire 
separation  of  the  flow  takes  place  when  AVEU  =  -.25.  This  is  not 
necessarily  so  in  the  general  case  of  an  axial-flow  compressor  stator 
channel  where  values  of  AVEU  <  -.25  have  been  obtained  without  flow 
separation.  It  would  be  expected  that  values  of  AVEU  <  -.25  could 
occur  far  downstream,  in  the  neighborhood  of  X  =  XMN,  and  that  calcula¬ 
tion  of  corner  flow  parameters  could  be  continued  from  this  point  on, 
assuming  no  flow  separation.  Likewise,  the  validity  of  the  results 
displayed  in  Figure  4  would  hold  and  the  limit  value  of  -.25  could  be 
used  in  estimating  the  fillet  radius  from  Figure  5. 

H12  and,  therefore,  n  are  calculated  for  the  main  flow  along  the 
pressure  and  suction  corners.  The  preceeding  comments  on  the  limit 
values  of  AVEU  apply  to  H12  since  it  is  also  a  separation  parameter. 
Boundary  layer  separation  occurs  on  an  airfoil  when  1.8  <  HI 2  <  2.2. 

For  values  of  H12  =  1.74  and  above,  Equation  5  yields  values  of  3  =  1 
and  3  <  1  when  6/eo  =  2.  Since  3  cannot  be  <  1 ,  Reference  1  ,  6/eo  must 
be  allowed  to  assume  values  >  2.  Figure  6  shows  the  variations  of  HI 2 
and  6/eo  along  the  pressure  corner  for  the  example  chosen  in  this  report. 

It  must  be  noted  that  the  variation  of  6/eo  is  to  be  interpreted 
as  an  average  value  of  this  parameter  between  the  entrance  plane  and  the 
value  of  X  at  which  the  maximum  value  of  H12  occurs.  In  Figure  6,  the 
peak  value  of  6/eo  is,  in  fact,  the  minimum  value  which  6/eo  may  assume 
for  the  associated  value  of  HI 2  (Reference  1).  Consideration  of  the 
experimental  values  obtained  for  3»  6/eo,  and  n  in  the  case  of  a  mildly 
adverse  pressure  gradient  (Reference  1)  makes  the  assumption  that  3 
is  very  near  1  for  the  strongly  adverse  pressure  gradients  prevailing 
in  compressor  stator  channels  quite  reasonable. 
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SECTION  III 
NUMERICAL  EXAMPLE 


The  numerical  data  used  in  this  example  was  obtained  from  Technical 
Report  AFAPL-TR-76-92  and  generated  with  a  single  stage  axial-flow  com¬ 
pressor  operating  at  100%  speed  near  peak  efficiency.  The  data  describes 
average  values  of  entrance  and  exit  flow  parameters  of  the  stator  channel 
and  its  physical  dimensions.  In  addition,  the  main  flow  velocity  profile 
along  the  suction  and  pressure  side  corners  was  scaled  from  data  calcu¬ 
lated  by  the  three-dimensional  computer  program  "NANCY  1,0"  of  Ai Research 
Manufacturing  Company  of  Arizona. 

Listed  below  are  the  average  values  of  the  velocity  at  entrance  and 
exit  planes  and  the  physical  dimensions  of  the  channel. 


Channel  width 
Channel  length 


2.54  cm. 
6.35  cm. 


Reynolds  number 
Meridional  Velocity 
Initial  Tangential  Velocity 
Exit  Tangential  Velocity 
Initial  Momentum  Thickness 
Initial  HI 2 


204.52  m/s 
253.18  m/s 
9.4  m/s 
.0038  cm. 


1  in. 

2.5  in. 

6 

10 

671 .  ft/sec 
831 .  ft/sec 
30.82  ft/sec 
0.0015  in. 

1 .47 


The  results  of  corner  flow  boundary  layer  parameters  calculation  as 

performed  by  the  computer  program  listed  in  the  Appendix  are  presented  in 

Figures  7  through  17.  From  these  results,  the  evaluation  of  the  corner 

fillet  radius  is  easily  performed.  For  the  pressure  corner,  for  instance, 

the  maximum  value  of  <5^,,  ,  which  occurs  at  S/L  =  .836  in.,  is  .037  (Figure 

-3/1-  6 

16).  It  is  seen  that,  from  Figure  4,  for  RE  =  10  ,  a  safe  value  of  63^ 
may  be  taken  as  0.007,  so  that  the  percent  decrease  in  sought  is 
81.2%.  Figure  5  is  then  used,  in  conjunction  with  the  average  Euler 


9 


AFWAL-TR-81 -2020 


number,  in  this  case  equal  to  -.246,  to  find  the  value  of  the  fillet  radius, 
which  is  0.22  in.  Similarly,  for  the  sunction  side.  Figure  17,  indicates 
that  the  maximum  value  of  <$3^  is  .04083,  occurring  at  the  trailing  edge. 

The  desired  decrease  in  63^  sought  is  82.86%,  which,  in  conjunction  with 

an  average  Euler  number  of  -.246  yields  a  fillet  radius  of  0.28  in. 

1 .  REMARKS 

The  discussion  presented  in  the  chapter  CONCLUSION  of  Reference  1 
applies  directly  to  the  contents  of  this  report.  This  is  particularly 
true  for  the  considerations  governing  the  fillet  radius  profile  along 
the  corner. 


10 


AFWAL-TR-81 -2020 


APPENDIX 

NUMERICAL  METHODS  (PROGRAM  FILLET) 

Tnis  program  involves  the  solution  of  three  linear,  first  order 
partial  differential  equations.  These  equations  (3,  4,  and  6)  are  inte¬ 
grated  using  a  fourth  order  Runge-Kutta  method.  The  parameter  $  of  the 
power  law  velocity  is  obtained  from  Equation  5,  using  a  polynomial 
approximation. 

1 .  PROGRAM  STRUCTURE 

This  program  consists  of  a  main  program  and  four  subroutines. 

(1 )  Main  Program: 

(a)  Calculates  the  free  stream  velocity  and  velocity  gradient 
along  a  streamline. 

(b)  Calculates  simultaneously  the  friction  coefficient,  endwall 
momentum  displacement  thickness  form  coefficient  HI 2,  and  corner  momentum 
displacement  thickness. 

(2)  Subroutines: 

(a)  Subroutines  "POLYNO"  and  "LINEQ"  use  a  polynomial  approxi¬ 
mation  to  calculate  the  free  stream  velocity  and  the  velocity  gradient 
along  the  suction  and  pressure  corners. 

(b)  Subroutine  "INTICK"  calculates  the  corner  interference 
displacement  thickness. 

(c)  Subroutine  "FILLET"  calculates  the  optimum  fillet  radius 
along  the  suction  and  pressure  corners. 

(d)  Subroutine  "AVEREUL"  calculates  the  coordinates  of  WSX 
minimum  used  in  the  calculation  of  the  Euler  number. 
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2.  PROGRAM  AND  SUBPROGRAM  DESCRIPTIONS 
Main  Program 

Calculation  Equation 

Determine  free  stream  velocity  1 

Determine  free  stream  velocity 

gradient  2 

Calculate  f°  (Ue/UQ)3'333  da 

Jo 

using  Simpson's  numerical 
integration  method 


First  Card 
103,  115 

111,  118 


149 


Determine  2-D  boundary  layer 
thickness  using  Truckenbrodt1 s 

quadrature  3  175 


Calculate  shape  factor  H12  4 

Calculate  corner  boundary  layer 
momentum  thickness  6 

Calculate  turbulent  3-D  endwall 
momentum  thickness  8 


198 

216 

220 


Subroutine  INTICK 

Calculation  Equation  First  Card 

Determine  interference  displacement 

thickness  9  5 

Subroutines  P0LYN0  and  LINEQ 

Calculation 

Calculate  M  +  1  coefficients  of 
polynomial  of  order  n  used  in  the 
calculation  of  WSX  and  DWSX. 


Subroutine  FILLET 
Calculation 


Calculate  optimum  corner  fillet 
radius. 
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Subroutine  AVEREUL 
Calculation 

Calculate  minimum  value  of  WSX  and  X 
at  which  this  minimum  occurs. 


3.  FORTRAN  NOMENCLATURE 


Variabl  e 
ALPHA 
AVEU 
BETA 

C,  CM 
CBL 

CFX 

CFX1 

CFX2 

CONSW 

DEL 

DELPRT 

DELX 

DEL2X 

DEL2XS 

DEL2X1 

DEL2X2 

DEL3 

DEL3MAX 


Descri ption 

Supplementary  angle  of  dihedral 

Average  value  of  Euler  number 

Parameter  governing  the  inclination  of 
isovels  to  the  wall 

Velocity  polynomial  coefficients 

Coordinate  of  intersection  of  bisector 
with  corner  boundary  layer  surface 

Skin  friction  coefficient 

CFX  at  X  +  DELX/ 2 

CFX  at  X  +  DELX 

X  component  of  free  stream  velocity  along 
mid-channel  streamline 

Extent  of  three-dimensional  corner  flow 

6/eo 

Incremental  step  along  X 

Boundary  layer  momentum  thickness  from 
Truckenbrodt1 s  quadrature 

Endwall  boundary  layer  momentum  thickness 
from  Suter's  hypothesis 

DEL2X  at>  X  +  DELX/2 

DEL2X  at  X  +  DELX 

Corner  interference  displacement  thickness 
Maximum  value  of  DEL3 
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Variabl e 

Description 

DEL3X 

Corner  boundary  layer  momentum  thickness 

DEWY 

Channel  entrance  and  exit  tangential 
velocity  difference 

DWSX,  DWSXM 

Free  stream  velocity  gradient 

EPSO 

Endwall  boundary  layer  thickness 

ETNA 

Distance  from  stator  mid-channel  endwall 
streamline 

EULE 

Euler  number 

EUM 

Average  value  of  Euler  number 

EXN 

Inverse  exponent  in  power  law  for  primary 
velocity  profile 

HI  2 

Form  factor 

NDA 

Corner  interference  displacement  area 

PERCNT 

Percent  decrease  in  DEL3MAX  required  to 
prevent  corner  flow  separation 

RAFIL 

Value  of  corner  fillet  radius  required  to 
obtain  percent  decrease  in  DEL3MAX 

SIBX 

Sine  of  angle  measuring  direction  of  free 
stream  velocity  W-R-T-  Y  axis 

SIGM 

X  values  for  VEL 

ST 

Length  of  channel  measured  along  mid 
stream! ine 

VEL 

Free  stream  velocity  input  data 

WSX 

Free  stream  velocity 

WSXO 

Total  velocity  at  channel  entrance 

WSX0 

WSX  initial  value 

WY1 

Tangential  velocity  at  channel  entrance 

X 

Distance  along  streamline 

YMIT,  YMIN 

Coordinates  of  WSX  minimum 
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4.  INPUT  DATA 

Three  input  data  formats  are  used:  alphanumeric,  integer,  and  real. 
Each  record,  or  card,  contains  only  one  type  of  data,  and  normal  FORTRAN 
convention  indicates  which  are  the  integer  quantities.  The  al phanumeric 
format  is  used  for  the  title  cards  and  consists  of  eighty  characters 
starting  in  the  first  column.  The  integers  are  entered  in  fields  of 
four  characters  with  up  to  twenty  numbers  per  card.  No  decimal  point 
may  be  used  and  the  number  should  be  right  justified  within  the  allocated 
field.  The  real  numbers  are  entered  in  fields  of  twelve,  with  up  to  six 
numbers  per  card  and  a  decimal  point  should  be  entered  to  ensure  correct 
interpretation  of  the  data.  The  first  field  on  each  card  begins  in 
col umn  1 . 

In  the  following  chart,  one  line  corresponds  to  one  card,  except 
where  it  has  been  necessary  to  use  two  lines  to  complete  the  description 
of  one  card. 

CARD  1 :  TITLE 

CARD  2:  X0,  XS,  ALPHA,  WSX0,  WSXO,  DEL2X1 

CARD  3:  DEL2X2,  DEL20,  HI  2 ,  ST,  WSFF 

CARD  4:  N0P,  N0D,  N,  M,  Ml,  L 

CARD  5:  WY1 ,  DEWY,  C0NSW,  ETNA,  DEL3X 

CARD  6:  VEL(K),  SIGM(K)}  This  card  occurs  M  times 

CARD  7:  VELM(K),  SIGM(K)}  This  card  occurs  Ml  times 

5.  DEFINITION  OF  INPUT  DATA  ITEMS 

TITLE  Alphanumeric  title  for  job.  Up  to  eighty 

characters  may  be  used  starting  in  column 

1. 

X0  Initial  value  of  left  hand  boundary  used 

in  subroutine  "AVREUL". 

XS  Initial  value  of  independent  variable 

"X"  along  a  given  streamline. 
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ALPHA 

WSX0 

wsxo 

DEL2X1 

DEL2X2 

DEL20 

HI  2 
ST 


Supplementary  angle  of  dihedral 
(Reference  1 ) . 

Initial  value  of  WSX. 

Total  velocity  at  channel  entrance. 

DEL2X  at  X  +  DEtX/2. 

DEL2X  at  X  +  DELX 

Initial  endwall  boundary  layer  momentum 
thickness  from  Suter's  hypothesis  (see 
note) . 

Form  factor. 

Length  of  channel  measured  along  mid 
streamline. 


WSFF  Channel  exit  free  stream  velocity. 

N0P  If  N0P  =  0,  Suter's  equation  1  is  used  to 

calculate  the  free  stream  velocities  along 
the  mid-channel  and  suction  or  pressure 
corner  streamlines,  using  the  following 
decimal  variables:  WY1 ,  DEWY,  C0NSW, 

ETNA. 

If  N0P  +  0,  a  polynomial,  1',  is  used  to 
calculate  the  non-dimensional  free  stream 
velocities  along  the  mid-channel  and 
suction  or  pressure  corner  streamlines 
from  the  experimental  data  entered  on 
input  data  cards  6  and  7. 

N0D  If  N0D  =  0,  all  computed  data  is  listed. 

If  N0D  f  0,  only  DEL3MAX,  EUM,  PERCNT, 
and  RAFIL  are  1 isted. 

N  Integer  variable  giving  iteration  interval. 

M  Number  of  velocity  data  points  along  corner 

stream! ine. 

Ml  Number  of  velocity  data  points  along  mid¬ 

channel  streamline. 


L  Degree  +  1  of  polynomial  used  in  the 

interpolation  of  velocity  data  points. 
Maximum  value  of  L  is  10. 
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WY1 

Tangential  velocity  at  channel  entrance 
(non-dimensional). 

DEWY 

Channel  entrance  and  exit  tangential 
velocity  difference  (non-dimensional). 

C0NSW 

X  component  of  free  stream  velocity  along 
mid-channel  streamline  (non-dimensional). 

ETNA 

Distance  from  stator  mid-channel  endwall 
streamline  (dimensionless). 

DEL3X 

Initial  corner  boundary  layer  momentum 
thickness  (see  note). 

VEL(K) 

Corner  free  stream  velocity  data  input. 
The  maximum  value  ;of  K  is  1 00 . 

SIGM(K) 

Distance  along  streamline. 

VELM(K) 

Mid-channel  free  stream  velocity. 

NOTE:  The  initial  value  of  the  two-dimensional  boundary  layer  momentum 
thickness  is  obtained  from  Falkner's  equation  (reference 
pp.  632) : 

cf  .0065 

2  (WSXO  *  DEL20/v)l/6 

with  the  friction  coefficient,  C^,,  obtained  .from  reference 
p.p.  600,  for  the  Reynolds  number  of  interest. 

6 

In  the  example  given  here,  Cf  =  .0045,  Re  =  10  ,  so  that  DEL20 
=  .00145.  T 

The  initial  value  of'  DEL3X  is  directly  obtained  from  the  initial 
value  of  DEL20,  using,  P.  Suter' s  assumption  that,  at  the  channel 
entrance,  the  ratio  DEL3X/DEL20  «  1.55 . 
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6.  OUTPUT  DATA 

The  printed  output  consists  of  two  sections.  First,  the  input  data 
is  completely  listed.  Following  this  are  printed  the  results  computed  by 
the  program.  If  N0D  =0,  the  bulk  of  the  calculated  data  is  printed  out. 

If  N0D  f  0,  only  selected  items  from  the  computed  data  are  printed 

out. 

Description  of  calculated  output  data: 

I.  N0D  =  0. 

Are  printed  in  the  following  order: 

The  coefficients  of  the  polynomial  used  in  the  interpolation 
of  WSX  data. 

The  calculated  values  of  WSX  at  data  points,  using  curve  fit 
polynomial . 

The  coefficients  of  the  polynomial  used  in  the  interpolation  of 
WSXM  from  data. 

The  calculated  values  of  WSXM  at  data  points,  using  curve  fit 
polynomial . 

The  minimum  value  of  WSXM  along  the  mid-channel  streamline. 

The  calculated  boundary  layer  parameters:  X,  WSX,  AVEU,  HI 2 , 
DEL2X,  DEL3X,  DEL2XS,  EPS0,  EXN,  BETA,  DEL3. 

The  required  fillet  radius  along  the  corner. 

The  percent  drop  in  DEL3MAX  to  prevent  separation. 

The  average  value  of  the  Euler  number  along  the  mid-channel 
streaml ine. 

The  maximum  non-dimensional  interference  displacement  thickness. 

II.  N0D  f  0 

Only  the  last  four  items  of  the  preceding  list  (N0D  =  0)  are 

printed. 
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Figure  1.  Non-Dimensional  Velocity  Along  Mid-Channel  and 
Pressure  Corner  Streamlines  Assuming  Circular 
Arc  Blade 
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Figure  3.  Influence  of  Endwall,  Boundary  Layer  Secondary 

Flow  on  the  Development  of  Corner  Boundary  Layer 
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Figure  4.  Interference  Displacement  Thickness  of  Turbulent 
Boundary  Layer  in  Convergent  or  Divergent  Flow 
Along  a  Corner  From  (3) 
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Figure  9.  Average  Euler  Number  Calculated  Along  Mid- 
Channel  Streamline 
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